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We study the P — >■ 77* (P = tt", j), r)') transition form factors by means of the local-duality (LD) 
version of QCD sum rules. For the case of ri and 7;', the conventional LD model provides a good 
description of the existing data. However, for the tt form factor we find disagreement with recent 
BaBar results for high even though the accuracy of the LD approximation is expected to increase 
with . It remains mysterious why the r) and r)' form factors to virtual photons, on the one hand, 
and the tt form factor, on the other hand, show a qualitatively different behaviour corresponding to 
a rising with violation of local duality in the pion case. In a quantum mechanical example we 
show that, for a bound-state size of about 1 fm, the LD sum rule provides an accurate prediction 
for the form factor for > a few GeV^. 

PACS numbers: 11.55.Hx, 12.38.Lg, 03.65. Ge, 14.40.Be 



1. INTRODUCTION 



The form factor describing the two-photon transition of a light pseudoscalar meson P is one of the simplest hadronic 
form factors in QCD. The corresponding amplitude 

{l{(li)l{(l2)\P{p)) = ie^^e^q^q.,_Fp^^^{ql,ql) (1.1) 

contains only one invariant form factor, Fp^^^{q\,q2). We shall be interested in the case of one real and one 
virtual photons, q\ = Q and — g| = > 0, and define Fp^{Q^) = Fp^^^{q\ = O,*/! = — Q^). For the pion case, the 
form factor FTr^{Q'^) has the following properties: (i) In the chiral limit of massless quarks and a massless pseudoscalar 
TT, the form factor at = is given by the axial anomaly [H, F^^{Q^ = 0) = 1/ {2^/2^:'^ f^), = 130 MeV. (ii) At 
large Q^, perturbative QCD (pQCD) predicts the asymptotic behaviour 0] F^^-yiQ^ 00) — V^f-n/Q'^- Brodsky and 
Lepage proposed a simple interpolating formula between these two values Q 

This formula may not provide an accurate description of the form factor at small nonzero Q but one may expect its 
accuracy to increase rather fast at values larger than a few GeV^. In contrast to these expectations, the BaBar 
collaboration recently reported a surprising result for the behaviour of the i^^^ form factor the product Q'^ Ft^^{Q'^) 
does not saturate at large but increases further (Fig. [T]). 
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Fig. 1: Form factor Ftt-/ as function of Q^: experimental data from ^ (red) and [1,0] (black); the solid line represents the 
interpolation (|1.2[) . 
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This result has aheady attracted a lot of attention in the literature (see, e.g. [6| and references therein). The aim 
of our analysis is to study the P7 transition form factor for P = 7r°, 77,77' by making use of the local-duality version 
of QCD sum rules . This approach allows one to study hadron form factors without knowing subtle details of their 
structure and to consider on equal footing form factors of different hadrons. 

We find that for both 77 and 77' transitions, the LD sum rule provides a satisfactory description of the form factors 
in the region ~ 5~ 100 GeV^ 0, The result from the LD approximation for the tt case, on the other hand, 
does not agree with the BaBar data. 

In order to test the accuracy of the LD sum rule for the P7 transition form factor, we consider, in parallel to QCD, 
a quantum-mechanical potential model. In the latter case, the form factor can be obtained both by the LD sum rule 
and by an exact calculation. Comparing these results with each other provides a probe of the LD approximation [lo| . 
Here we find that, independently of the details of the confining potential, the LD sum rule reaches the accuracy of a 
few percent already at relatively low values of the momentum transfer. 

It remains mysterious why, in contrast to the success of the LD approximation for the 77 and 77' transitions and to the 
experience from quantum mechanics, the results from the LD approximation strongly contradict the BaBar result. 

This Letter is organized as follows: In the next section, we recall the structure of the (VAV) amplitude and the 
relation of F^^j to the axial anomaly. Section 3 presents the dispersion representation for (VAV) and touches the 
issue of multiloop radiative corrections to this quantity. Section 4 reviews the LD approximation for the P7 form 
factor and studies the expected accuracy of this approximation making use of a quantum-mechanical testing ground. 
Section 5 applies the LD sum rule to the Fp-y, P = 7r°,?7,7;' form factors. Section 6 summarizes our conclusions. 



2. THE THREE-POINT FUNCTION (VAV) AND THE AXIAL ANOMALY 



Let us start with the amplitude of two-photon production from the vacuum induced by the axial-vector current of 
arly massless quarks of one flavour, = q^^'y^q, with £1^2 denoting the photon polarization vectors: 

(7(9i)7(92)|j'(x = 0)|0) =T^c«/3(pkl,g2)£?£^, P = qi + q2. (2.1) 



The amplitude is obtained from the vacuum expectation value of the T-product of two vector and one axial- 

vector currents and will be referred to as the (VAV) amplitude. Vector-current conservation yields the following 
relations: 

T^,Mp\quq2)q? = 0, T^MpIQu Q2)ql = 0. (2.2) 
The general decomposition of the amplitude contains four invariant form factors and may be written as 

Tfj,al3{p\qi,q2) = n 2 ^al3qiq2^FA + {ql£tj.aflq2 — qiait^qil3q2 2 2 ^1 '=9lQ;3g2 

p^ — nip f"^ _ j^i^^ 

+ {q2^l-il3aqi - q2l3£fj.q2aqi 2^^ 2 qhq2l3cyqi)iP2 + (qi - 92 )^iea/3<;ig2 ^-f's ■ (2.3) 

p ni p 

The explicit calculation of these form factors at one and two loops in QCD leads to F3 = [llj. Therefore, we 
shall omit the corresponding term. The absence of any contact terms in T^a^ can be verified by reducing out one of 
the photons and using the conservation of the electromagnetic current . 

The parameterization (12. 3p takes into account the pole at p^ = TTip, related to the contribution of the lightest 
pseudoscalar state. Apart from the pole at p^ = Trip, the amplitude has no singularities at small p^. In the chiral 
limit, TTip = 0, both the second and third Lorentz structures in (|2.3p are transverse with respect to p^, and the 
form factor Fa represents the axial anomaly [l| . According to the Adler-Bardeen theorem [l^ , the axial anomaly is 
saturated by the one- loop expression, Fa = l/(27r^), and remains non-renormalized by higher-order corrections. 

Separating the longitudinal and the transverse structures for the case TTip 7^ 0, we obtain 

Tf,afl{p\qi,q2) = -^eafiqiq2iFA 2\^P^»l3qiq2'iiFA + qjFi + qlF2) 

P P \P '"T^pl 

I / 2 Pt^ 2 \- IT 

+ [qi<^fial3q2 - qiaf^tiqil3q2 " 9l ^91 a/3g2 1 

+ (92£/i/3a(;i - q2p£p,q2aqi - ■^'72<^92/3Qgi )*^2 • (2.4) 

Beyond the chiral limit, the first two terms contain two poles [l3|: a "kinematical" pole at = 0, which cancels the 
corresponding singularities in the transverse Lorentz structures, and the "dynamical" pole at p^ = m?p corresponding 
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to the TT meson. The full amplitude is regular at = 0. That is, the pole at = in the transverse Lorentz structures 
is of purely kinematic origin and does not correspond to a massless particle. By forming the divergence, one obtains 



92 



Fa ~ 2 ^ {Fa + qiFi + qiF2) 
nip — 



(2.5) 



As is clear from (j2.5p , the transition form factor of interest is given by the linear combination of Fa , Fi , and F2 



fl 



-{Fa 



qlFi 



qlF2)\p2 



(2.6) 



Thus, the P — > 77 form factor is proportional to the axial anomaly only at one kinematical point, qf ^ q2 = 0. 
At this point, the pion form factor in the chiral limit is protected from radiative corrections by the Adler-Bardeen 
theorem; the one-loop result represents the exact result. Hence, in the chiral limit one expects the pion pole at = 
to emerge in one-loop diagrams for the (VAV) amplitude. Indeed, in this unique situation a pole dual to a single 
hadron state emerges from the single one-loop diagram of perturbation theory |15| . 

However, already if one of the photons is virtual, the transition form factor Fp^^^ and the axial anomaly are not 
proportional to each other. It is the topic of the next two sections to analyze the form factor Fp^{Q'^) by means of 
dispersive sum rules [l6l |. 



3. DISPERSION REPRESENTATIONS FOR (VAV) AND THE AXIAL ANOMALY 



We now discuss the one-loop expression for the amplitude. To this end, a slightly different parameterization, 
obtained by setting {Fa + qfFi + (Z2-^2)/(p^ — rrip) — Fo, proves to be convenient: 



/jgJiPl -I- {qj )iF2 



As follows from ()2.6|) . Fq contains the contribution of the pseudoscalar meson of our interest. 
We also consider the transition amplitude of the pseudoscalar current operator 9759: 

(7('7i)7(92)|g75'7|0) = eai3qiq2£i4F5{ql,ql,p'^)- 
The two-photon amplitude of the divergence of the axial current takes the form 

{l{qih{q2Wjl\0) = e^pg,g,e^el{p^Fo - qjF, - 91^2). 



(3.1) 



(3.2) 



(3.3) 



The case of our interest is qf = 0, then the form factor Fi does not contribute to the divergence. In perturbation theory. 



the form factors Fq^ F2, and F5 may be written in terms of their spectral representations in p (with q 
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F^{p'.q') 



■Ms,q')- 



(3.4) 



4m2 



To one-loop order, the spectral densities read [12, llTHli 

1 1 



Ao(s,g2) 
A2(s,g2) 
A5(s,g2) 



2n {s 
1 



1 



2n (s-g2)2 
1 m 



2tt s — q- 



■log 



-q w + 2m log 

+ 2m^ log 
1 + w 



1 + w 
1 ~ w 

1 -I- u; 

I — w 



1 — w 



w = a/1 — 4m2/s. 



(3.5) 



Obviously, the absorptive parts obey the classical equation of motion for the divergence of the axial current 

sAo{s,q^) -q^ A2{s,q^) = 2mA^{s,q^). (3.6) 

The form factors then satisfy 



p'Fo{p',q')-q'F2{p',q') - 2mF^{p\q') 



ds Ao(s, q^). 



(3.7) 
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The last integral is equal to — l/27r, independently of the values of m and g^, and represents the axial anomaly 

p^Foip^q^) - q'F2{p^q^) = 2m F^ip^, q') + (3.8) 

In the chiral limit m = and for q^ = 0, the form factor Fq develops a pole related to a massless pseudoscalar meson 
[isj . The residue of this pole is again the axial- anomaly l/27r^. 

As is clear from p.7p . the anomaly represents the integral of Aq, the spectral density of the form factor Fq. Adler 
and Bardeen tell us that the anomaly is non-renormalized by multiloop corrections. The easiest realization of this 
property would have been just the vanishing of multiloop contributions to the spectral density Ao(s, q"^). An argument 
in favour of this possibility comes from explicit two-loop calculations [ill , fisj ] which report the non-renormalizability 
of the full {VAV) vertex to the two-loop accuracy. However, if so, the full form factor Fq is given by its one-loop 
expression. Then, this expression should develop the pion pole, known to be present in the full amplitude for any 
value of q^. But, obviously, this pole does not emerge in the one- loop expression for i^o if 9^ 7^ 0! 

This requires that multiloop corrections to the form factor Fq (and, respectively, to its absorptive part Ao(s,q^)) 
do not vanish. Then one may ask oneself how it may happen that the anomaly nevertheless remains non-renormalized 
by multiloop corrections? The only possible answer we see is that the non-renormalization of the anomaly is reached 
due to some conspiral property of multiloop contributions to Ao(s,q^) forcing their integral to vanish (in fact, quite 
similar to the one-loop result: although the spectral density explicitly depends on m and q^, its integral is an m- and 
g^-independent constant). 

So we conclude that, in spite of the fact that explicit calculations yield a non-renormalizability of the full (VAV) 
vertex to two- loop accuracy [HI, [ll], the conjecture of PJJ that this result might hold to all orders of the perturbative 
expansion may not be valid. 



4. Fp^ FROM A LOCAL-DUALITY SUM RULE FOR Fo 

As is obvious from (|2.4p . the contribution of the light pseudoscalar constitutes a part of the form factor Fq. The 
Borel sum rule for the corresponding Lorentz structure reads (Q^ = —q^ > 0) 



/ 



(is exp(— sr)Ao(s, Q^) = —.fpFpj{Q^) exp(— 77ipT) + contributions of excited states. (4-1) 



Exploiting the concept of duality, the contribution of the excited states is assumed to be dual to the high-energy region 
of the diagrams of perturbation theory above an effective threshold Scff- After that, setting the Borel parameter t = 
(which yields the so-called local-duality limit), we arrive at the LD sum rule for a pseudoscalar ^g-meson 

dsAQis,Q^)^~fpFp^iQ^). (4.2) 

The spectral density Aq to one-loop order is given by p.Sp ; two-loop corrections were found to be absent [ll|, [ll| . As 
discussed above, higher-loop corrections to Aq cannot vanish; so the l.h.s. of (|4.2p is known to 0{al) accuracy. Recall 
that all details of the nonperturbative dynamics are encoded in a single quantity, the effective threshold Scff(Q^)- The 
effective threshold is an essential parameter of the method of dispersive sum rules; it is not identical to the physical 
threshold and is not universal (i.e., it is specific for the correlator under consideration); moreover, in general it depends 
on Q2 ^20]. 

In the chiral limit, the LD expression for the form factor for the one-flavour case is particularly simple: 

Fp.iQ')-A-r JF^^- (4-3) 

Apart from neglecting and higher-order corrections to the spectral density Aq, no approximations have been done 
up to now: we have just considered the LD limit r = 0; for an appropriate choice of s^siQ"^) the form factor may still 
be calculated exactly. Approximations come into the game when we consider a model for Soff(Q^). 

Irrespective of the behaviour of Soff(Q^), at = the form factor is related to the axial anomaly: Fp^{Q) = 
l/(27r^/p). QCD factorization requires ScsiQ"^) — > 47r^/|, for large . The simplest model compatible with this 
requirement is obtained by setting 

Scff(g')=4^Vp (4.4) 
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for all values of [2l[. We shall refer to the choice (I4.4p yielding for the neutral pion case the Brodsky-Lepage 
result (|1.2p . as the "conventional LD model". ^ The relevance and the expected accuracy of the LD model may be 
tested in those cases where the form factor Fp^{Q'^) is known, i.e., may be calculated by other theoretical approaches 
or measured experimentally. Then, the exact effective threshold may be reconstructed from (|4.3p . in this way probing 
the accuracy of the LD model. Let us play this game making use of a quantum-mechanical model, where the exact 
form factor may be calculated from the solution of the Schrodinger equation, and then compared with the result of 
the sum rule for a three-point function of nonrelativistic field theory. 



A. LD sum rule in quantum mechanics 

The analogue of the 7r7 form factor in quantum mechanics is given by (23l . [23 | 

oo 

FnK{q) = J dT{^\J{q)exp{^HT)\r = 0), (4.5) 



where the current operator J{q) is introduced by the kernel {r'\J{q)\r) = exp{iq ■ r) S^-^\r — r'), the Hamiltonian H 
governing the nonrelativistic potential model reads 

H=^ + V{r), (4.6) 

and ^I' is the corresponding ground state. Technically, the LD model for the form factor (|4.5p is constructed from the 
quantum-mechanical analogue of the three-point function in the same way as for the case of the elastic form factor (for 
details, consult [id])- Recall, however, an essential conceptual difference between the form factor and the elastic 
form factor with respect to the factorization of these quantities at large momentum transfers: The factorization of the 
elastic form factor requires the presence of both Coulomb and confining terms in the interaction. The factorization 
of the Pj form factor does not require the presence of a Coulombic term and emerges also for a purely confining 
interaction. The LD model for a given form factor is tightly related to its factorization properties; specifically, the 
LD model for the analogue of the P7 form factor may be formulated in quantum mechanics for the case of a purely 
confining potential. 

Figure [2] depicts the exact effective threshold kcs{Q) — the quantum-mechanical counterpart of the effective thresh- 
old Seff (Q^) — for the example of the harmonic-oscillator potential V{r) = moj'^r'^ /2, for parameter values relevant 
for hadron physics: a reduced mass of the light quark of m = 0.175 GeV and an interaction strength of w = 0.5 
GeV, which lead to a size of the ground state around 1 fm, a typical size of a ground-state hadron in QCD. We have 
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Fig. 2: Effective threshold kcsiQ) for the three-point function in quantum mechanics, recalculated from the exact form factor 
-Fnr(Q), Eq. (|4.5p . in the harmonic-oscillator potential model. Rg = |^'(r = 0)^. 



^ In an alternative approach to the P7 form factor i22j] . the pseudoscalar meson is described by a set of distribution ampUtudes of 
increasing twist which are treated as nonperturbative inputs. In our analysis, the deviation of the effective threshold Sofr(Q^) from its 
asymptotic value An'^fp corresponds to some extent to the contribution of higher- twist distribution amplitudes in the approach of |22|| . 
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checked that a similar picture for the effective threshold emerges for other confining potentials; moreover, adding the 
Coulomb potential changes this picture only slightly. 

From the behaviour of kcs{Q), we conclude that the LD model may be expected to work increasingly well already for 
above a few GeV^. Inspired by this result, we now look what the LD model predicts for (tt, 77, rj') — > 77* transitions. 

5. THE (7r°,T7,T7') ^77* FORM FACTOR 

A. The tt" — >^ 77* form factor 

Taking into account the 7r° flavour structure and choosing the relevant interpolating current = {u^^lI^u — 
<^7/i75'^)/'\/2, the LD sum rule reads 



V2V9 9; 27r2/pSeff(Q2)+Q2' 



Nr. 



(5.1) 



which by setting Soff(Q^) = 47r^/^, /jr = 130 MeV, leads to the Brodsky-Lepage formula (11.21) . Figure [T] shows the 
corresponding plot. Figure |3] represents the "experimental" effective threshold recalculated from the form factor data 
via Eq. (|5.ip . This "experimental" effective threshold may be well approximated by a linearly rising function of 
. Surprisingly, the BaBar data show very strong — and growing with — violations of local duality! This 
observation is in absolute contradiction to our experience from quantum mechanics. Let us investigate next what 
happens in the case of the 77 and 77' mesons. 



Seff(Q^) [Gel/] 




10 15 20 25 30 35 40 



[Ge\/] 



Fig. 3: Effective threshold Soff((5^) recalculated from the data Q (red) and 0, (black) by means of the LD relation (|5.ip for 
the form factor Fn-y{Q^). 



B. The (ri,ri') — >^ 77* form factor 

The simple expression (|1.2I) is sometimes erroneously assumed also for the 77 and 77' cases [1, H^. However, the naive 
replacement /^r — > /t),j)' in (II. 2p yields wrong expressions for Fi-^^/)^. The correct way to proceed is to take into account 
the presence of two — nonstrange and strange — components in the 77 and rj' mesons and their mixing. Making use of 
the 77-77' mixing scheme from [2g, [27[ (see also [l^), the flavour structure of 7/ and rj' may be described as follows 



,uu + dd. 



\ss) snif 



V2 



I smffl + \ss) sin( 



The corresponding expression for the form factors take the form 

V2\9 9 



Fri'jiQ ) 



F„(Q2)sin( 



39.3", 



^F,{Q^) sinq 
^F,(Q2)cos< 



(5.2) 



(5.3) 
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Here Fn{Q'^) and Fs{Q'^) are the form factors describing the transition of the nonstrange and ss-components, respec- 
tively. The corresponding LD sum-rule for these form factors have a simple form 



Jn J 


Fs^iQ'') = J j dsA,{s,Q^). (5.4) 



A„ and correspond to Aq with different quark masses in the loop. In numerical calculations we set m,i = rud ~ 
and nis = 100 MeV. Accordingly, the LD model involves two separate effective thresholds for the nonstrange and the 
strange components [27l [: 

=47r2/2, /„«1.07/,, si'^ =^7T'fl /.«1.36/,. (5.5) 

The LD model may not perform well for small values of , where the true effective threshold is smaller than the LD 
threshold. However, for larger the LD model gives reasonable predictions for the form factors, as illustrated 
by Fig. H 



Q^F„(Q') [GeV] 



Q^F^.^iO') [GeV] 



0.25 










0.2 












0.15 












0.1 
0.05 


f ^ 





5 10 15 20 25 30 35 40 



lGe\/] 




♦ + 



TT 



10 15 20 25 30 35 



[Gel/] 



Fig. 4: LD form factors Fr/y and Fn'-y vs. Q^, compared with the experimental data presented in [1,01 (coloured dots) and 
[§| (black dots). The data points "borrowed" from the timelike momentum transfer = — — 112 GeV^ [1], q'^Frjj{q'^) — 
0.229 ± 0.03 ± 0.008 GeV and q^F^-^i{q^) = 0.251 ± 0.019 ± 0.008 GeV are not shown in this plot. 



6. CONCLUSIONS 



In this paper, we analyzed the P — >■ 77* transitions for P ~ tt^ , rj, rj' . 

I. We emphasized that the P — 77* form factor is proportional to the axial anomaly only if both photons are on-shell; 
if at least one of these photons is virtual, this proportionality is lost. As a result, the P — >■ 77* form factor in the chiral 
limit is not protected from receiving higher-order radiative corrections by the Adler-Bardeen theorem. Moreover, for 
virtual photons, the one-loop expression for the (VAV) amplitude does not develop a pole at = 0, related to a 
massless pseudoscalar. Therefore, the one-loop result for the form factor Fq cannot represent the full result for this 
quantity. Fq should receive radiative corrections at higher orders in the loop expansion, in spite of the absence of 
two-loop radiative corrections to the (VAV) amplitude reported in [111, [13 ■ 

Interestingly, the form factor Fq and the axial anomaly are given by dispersive integrals involving the same function, 
Ao(s,g^). By the argument given above, radiative corrections (coming from three and more loops) to the spectral 
density Ao(s,(3'^) cannot vanish. Then, the Adler-Bardeen theorem requires some specific conspiral properties of 
multiloop corrections to Ao(s,(j^) enforcing the vanishing of their integrals over s. 

II. We applied the local-duality version of QCD sum rules to the {7T°,ri,r]') — >■ 77* transition form factors. An 
attractive feature of this approach is the possibility to study form factors of bound states without knowing subtle 
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details of their structure. Moreover, it allows one to consider on equal footing form factors of different bound states. 
Our findings may be summarized as follows: 

1. We tested the accuracy of the LD model for the Pj transition form factor in quantum mechanics. We calculated 
this form factor from the solution of the Schrodinger equation and compared with the result of the quantum- mechanical 
LD sum rule. This comparison reveals that for a usual bound state, with a typical hadronic extension of about 1 fm, the 
LD sum rule is expected to yield accurate predictions for the form factor for larger than a few GeV^ ; this accuracy 
increases with Q^. At small but nonzero Q^, deviations from the LD model depend on subtle details of the confining 
interaction. 

2. Surprisingly, the BaBar data for the pion form factor exhibit an extreme violation of local duality in the 7r7 form 
factor even at = 40 GeV^. Moreover, the violation of local duality increases with in the range = 10—40 GeV^. 

3. Even more surprisingly — taking into account the strong disagreement in the pion case — the LD predictions 
agree with the experimental data for both the rj and rj' mesons in the rather broad range = A — 100 GeV^. 

The question why the nonstrange component in rj and 77', on the one hand, and the pion, on the other hand, should 
lead to a qualitatively different behaviour of the P — > 7 7* form factor remains mysterious. So far no compelling 
theoretical explanation for this strange phenomenon has been found (see also p2l. [29|). 
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